Surface electron mobility over a helium film 
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The mobility of surface electrons localized over helium films underlying solid substrates has been 
evaluated by solving the Boltzmann equation in the time relaxation approximation and the force 
balance equation in which an effective mobility is obtained in terms of the dynamical structure 
factor of the nondegenerate electron liquid. The essential processes of electron scattering by gas 
atoms, ripplons, and film-solid interface roughness are taken into account. The electron mobility 
dependence on the film thickness and temperature is determined and compared with experimental 
data available. We find that the interface-roughness scattering is the dominant process for explaining 
the experimental results. We estimate the extended defect sizes of the underlying substrate within 
the Gaussian correlated model for interface roughness. 

PACS numbers: 73.20.-5, 7323.-b,68.90.+g 



I. INTRODUCTION 

Surface electrons (SE) levitating over liquid helium has 
been used as a paradigmatic quasi-two dimensional elec- 
tron system (Q2DES). It constitutes the counterpart of 
the electron system on semiconductor heterostructures 
in the sense that the electrons on helium obey the clas- 
sical statistics because the low densities achievable in 
experiments (< 10 9 cm -2 ). One point difference be- 
tween them comes from the scattering mechanisms. In 
the former system, electrons are scattered by surface ex- 
citations, the ripplons, and by vapor atoms of the he- 
lium surface. In the latter one, scattering by impurities, 
phonons, and interface roughness are the processes that 
limit the electron mobility. Besides the remarkable phe- 
nomena discovered in this low-dimensional electron sys- 
tem, for instance, Wigner crystallization! and the exis- 
tence of edgemagnetoplasmons 2 ^ SE have been used as a 
probe to investigate the elementary surface excitations of 
cryogenic liquids and solids. Very recently there is an in- 
tensive search for experimental realization of SE as qubits 
leading to a quantum computer For a overview of the 
field see Ref. 

For electrons trapped on the bulk helium surface it is 
well known that the SE mobility is dominated by the 
electron-helium gas atom scattering for T > 1 K whereas 
the SE scattering by quantized capillary waves (ripplons) 
is responsible for mobility at lower temperatures where 
the helium vapor density becomes negligible. The ex- 
perimental data for SE mobility over bulk helium are in 
reasonable agreement with theoretical calculations. 8 

The situation changes when SE are floating over a he- 
lium film which in turn is deposited over a solid substrate. 
In such a condition it is possible to increase significantly 
the accessible range of electron densities- and the elec- 
tron correlations may become important. Furthermore, 
the electron energy spectrum, the ripplon dispersion rela- 
tion and the electron-ripplon interaction are modified sig- 



nificantly due to film effects and the van der Waals forces 
from the substrate now play a decisive role in the trans- 
port properties. Furthermore, besides the usual scatter- 
ing mechanisms pointed out above, we must consider the 
SE scattering by interface defects at the helium film - 
substrate boundary, which, as it have been shown in ex- 
periments is responsible by the unusual behavior of the 
SE transport properties on a helium fil m ,i2iiiii2ii&±i 

The present work intends to provide a detailed the- 
oretical study of the mobility of SE localized over a he- 
lium film deposited on substrates, as solid neon, glass and 
poly(methyl-methacrylate) [pmma], which are the mate- 
rials that have been used in the experiments. We em- 
ploy the Boltzmann transport equation approach (BEA) 
and the force balance equation method (FBEM) in which 
Coulomb effects on the mobility are taken into account 
through the dynamical structural factor of the Q2DES. 
We calculate the SE mobility as a function of tempera- 
ture and the film thickness in the range 100 < d < 1000 
A. From the comparison with experimental results we 
check the role of all scattering processes involved and ver- 
ify the reliability of the interface roughness model used in 
our treatment. We find that the SE mobility is strongly 
determined by the interface roughness substrate and the 
results are in quite good agreement with the experiments. 



II. BASIC RELATIONS 

A. Surface electron states 

Electron states on a liquid helium film are confined in 
the direction normal to the surface (z) due to the poten- 
tial well created by a infinitely high barrier which pre- 
vents the electrons penetrate inside the liquid phase and 
by and attractive potential due to the electron interaction 
with the polarizable substrates and an applied holding 
electric field E± along the z direction. The electron po- 
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tential energy over flat helium film located at — d < z < 
is well known and given byi^ 



z + d 



eE± 



(1) 



where A = e 2 (e He - l)/4(e ffe + 1) and Ai = e 2 £ He (e s - 
£He)/ (£He + l) 2 (e s + £He) with £ He and e s being the 
dielectric constants of helium and solid substrate, respec- 
tively, and e the electron charge. 

For electrons moving freely along the helium surface 
with the wave function and the energy spectrum given 
respectively by tfji(r,z) = A~ x l 2 exp(ik • r)xi(z) and 
Ei(k) = H 2 k 2 /2m + A/, where k and r are the wave 
vector and position vector in the plane (cc, y) of the liq- 
uid interface and A is surface area occupied by electrons. 
Unfortunately an analytical solution of Schrodinger equa- 
tion for Xi( z ) an d subbands energies A; cannot be found 
for U(z) given by Eq. Q). In this work, we use the 
variational method by choosing the trial wave function 
Xi(z) — 27 3 / 2 ze~ 7Z . The parameter 7 depends strongly 
on d and is determined from the subband energy mini- 
mization as was done in Refs. [Hand G3 One has found 
that the energy gap between the ground and first excited 
subbands increases considerably by decreasing d. For ex- 
ample A 2 — Ai ~ 12.7 K for d = 1CP 6 cm for the neon 
substrate and is much larger for a substrate with higher 
£ S *L Then one can disregard the possibility of electron 
transition to excited subbands which is proportional to 
exp[— (A2 — Ai)/T] and take only the ground subband 
into account in the calculation of the scattering matrix 
elements. 



B. Interaction Hamiltonians 

The electron interaction with helium vapor atoms by 
a contact type Hamiltonian is given by 
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(2) 



mind, we write the Hamiltonian for the electron-ripplon 
and electron-interface couplings in a unique form as 



(ei) 



(3) 



within 



V m {z) = — 



--Ktiqz) 
qz 



Ai 



{z + dy 



eEi 



and2a 



V lq {z) 



-Ai 



2e He \ qK ± [q (z + d)] 



+ 



e s + E H e J {z + d) 
£s - £He \ q 2 K 2 [q{z + d)] 
£ s + £Hc 



where Kj{x) is the modified Bessel function. Here we 
used the Fourier transform £j (r) = A~ x l 2 £jq eiq ' r 
for £1 (r) and £2 (r) and, quantizing the oscillations of 
free helium surface, one obtains 



' hq t&nh(qd) 



2pu q 



(4) 



with a q and a q being the annihilation and creation op- 
erators for ripplons with wave number q and satisfying 
the dispersion law 



(5) 



4 = ( -q 3 + g'q) tanh?rf, 



where a is the surface tension coefficient and p is the 
density of helium. In Eq. (JSJ, g' = g + 3(3/ pd 4 , with g 
being the gravity acceleration and (3 the van der Waals 
constant. 

Even though one cannot find from first principles the 
interface roughness displacement £2 (r), we will use a rea- 
sonable model to calculate the interface scattering con- 
tribution to SE transport properties. 



where R e and R a are three-dimensional positions of elec- 
trons and atoms, respectively and ao ~ 0.61 A is the 
scattering length^ The electron-ripplon interaction was 
derived in Ref. and the electron interaction with sur- 
face roughness was obtained in the case of SE over solid 
hydrogen. 19 

Very recently two of us have constructed an Hamilto- 
nian where both electron-ripplon and electron-solid in- 
terface interactions are treated on the same footingiSi 
The resulting total interaction potential for any substrate 
was obtained from the solution of the Poisson's equa- 
tion through the perturbation of the boundary condi- 
tions from the flat positions at z — and z — —d. Due 
to very small polarizability of liquid helium, the result 
of our approach for the electron-ripplon interaction co- 
incides with that obtained in Ref. Il5l Bearing this in 



III. TRANSPORT PROPERTIES 

A. Boltzmann kinetic approach 

In order to calculate SE mobility along the helium sur- 
face in the presence of a driving electric field E|| we ap- 
ply the well-known kinetic formalism where the electron 
scattering plays a crucial role in transport properties and 
does contribute to collision integrals in the Boltzmann 
equation as follows 



= S eg {/} + S er {/} + S ei {/} , 



(6) 



where /(k) is the electron distribution function and 
&eg {/}> S er {/}, and S e i {/} are collision integrals of 
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electron with helium atoms, ripplons and roughness in- 
terface between the film and the substrate, respectively. 
By considering quasielastic scattering processes the solu- 
tion of Eq. (7) is /(k) ~ / (0) (fc) + cos <p where 
f(°'(k) is equilibrium distribution given by the Boltz- 
mann function for the nondegenerate Q2DES, cp is the 
angle between k and E||, and 



' hv(k) dk ' 



where the collision frequency v(k) = v eg (k) + v er (k) 
u e i(k) and 



1 - 



k'.k 



(7) 



The probability amplitude per unit-time Wj(k, k') for 
electron transition from states k to k' is given by Fermi's 
golden rule and the subscript j corresponds to each scat- 
tering mechanism whose interaction potential is given by 
Eqs. and ©. The SE mobility is then given by 



e 
m 



xe x dx 
i (krx 1 / 2 ) 



(8) 



where x — k 2 /kj,, with the thermal wave vector kx — 

The frequency of electron-atom collisions has been cal- 
culated in Ref. 18 and the result is 



3n aghugj 
8m ' 



(9) 



Note that v eg does not depend on k being a func- 
tion of volume concentration of helium atoms n g — 
(MT/2Trh 2 ) 3 / 2 exp(-Q/T), where Q is the vaporization 
energy and M is the 4 He mass. As n g is an increasingly 
exponential function of T, v eg becomes negligible in com- 
parison with v er and v e d for T < 1 K. 

The calculation of v er and v e i is more complicated. Us- 
ing Eqs. an< i one can obtain in a straightforward 
way that both ^ er (fc) and v e i(k) can be written as 



2A- 



Ver(i) (k) 



m f q < |£l(2)q| > \{Xl\Vr(l)q{z)\Xl)\ 



rh 3 k 2 



coth(hu> q /2T) ~ TFjhwq ^> 1 for long wavelength rip- 
plons with huj q <C T. Equation H1(J|) results into 



v er {k T x 1 ' 2 ) = v^)J{x) 
where v^) = 870 T 2 /it ha and 

J(x) = 



(12) 




Here 



18 AnT 2 



$ er (y) = (j>(y) + 



Ai 



[1 + 2jd+ 



Aoy 2 

47d(l + 7rf)e 27d £;i(-27d)] , 



(14) 



with 

0(2/) 



1-y 2 



In 



vV 



9(1 -y) 

e(y-i), 



In the above expressions A = fi .7 /2m, Aj_ = 
h 2 -/l/2m,, A = mAg/2fi 2 , 7 _l = (3me£: ± /2fi 2 ) 1 / 3 , x c = 
H pg' /8maT, Ei(y) is the exponential- integral function, 
and 8(y) is the step function. 

For the scattering by interface roughness, one obtains 



16mA 2 fc 4 
irh 3 



Vei (k) 



v/4fc 2 - g 2 



-d<j, 



x$ 2 



7T/2 

y sin 6 6» (|6 9 (g = 2fcsin^ 



■sin^ dO 



(15) 



(10) 

where (...) means an average over ensemble. Considering 
the electron-ripplon scattering one obtains, from Eq. 



where 

(y) 



1 f x 2 Kx [y (x + 2jd)} 



x + 2jd 



exp(— x)g?x. (16) 



l(2)q| 



hq tanh(qe?) 
2puJ q 



{2N q + 1) 



(11) 



where the ripplon number N q is given by the Bose- 
Einstein function, which leads to 2N„ + 1 = 



In Eq. 116|) and thereafter Kj (x) is the modified j-order 
Bessel function. Now we must to choose an interface 
roughness model to derive the amplitudes £29 appearing 
in Eq. (|l(Jfl . Following Prange and NeeSI, we use the 
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Gaussian correlated model, in which the interface rough- 
ness is described by two characteristic sizes, the height £ 
and the lateral length I. These parameters define the au- 
tocorrelation function for interface roughness which can 
be written as 



(£(r)£(r')> = ^exp 



n'|2 



which leads to 



(17) 



(18) 



Assuming this model, the collision frequency, given by 
Eq. Ijl5(l . can be calculated and the result is 



(krx 1 ' 2 ) 



v<$xF 



(k^l) , (19) 



where vf) = 32m 2 A? $T/h 5 and 

tt/2 

F(b) = b 2 J sin 6 6>$ 2 4 f A sin oj exp (~b 2 sin 2 0) d6. 
o 

(20) 

The Gaussian correlated model was employed due to 
the successful description of the interface roughness for 
similar systems, for its simplicity and, since there are 
only two fitting parameters I and £o- 



B. Many-electron effects 

Monarkha and co-workers 2 ^* 2 ^ 2 ^ 2 ^ have studied elec- 
tron correlation effects in the transport as well quantum 
magnetotransport of SE on helium by using the force 
balance transport equation method (FBEM). In this ap- 
proach the frictional force experienced by the center of 
mass of the Q2DES due to electron-scatterer interactions 
is evaluated through the calculation of the momentum 
rate absorbed by the scatterers. For the system of N e 
electrons the kinetic frictional force F f r is given by 



7r 



= iV e eE|| = - 



dP 

dt' 



(21) 



where 

dP 

~dt 



^2 (P»' -Vv)\( v 'ij'\Hint\V,i) 
\v',j' 

8[{E U , + E r ) - {E„ + E,)]) . 



(22) 



Here p„ and E v are momentum and energy of scatterer 
system in the state \v), Ej the energy of the electron 
system in the state \j), Hi nt is the Hamiltonian of the 
interaction between the electron and the particular scat- 
terer. The angle brackets means a thermodynamic av- 
erage in the laboratory coordinate system. In Ref. 03 



is was shown that dP / dt can be written in terms of the 
dynamic form factor of the Q2DES 



s lab (q,c) = ^ ( l(i>-qli)l 2 <5 ( E o' - E i - M 



in the laboratory frame. <Siab(q, w) is connected with 
the form factor 5(q, u>) in the center-of-mass of the elec- 
tron system, moving at a velocity u relative to the lab- 
oratory frame, as 5i a b(q, w) — S(q,u — qu). Here 
ri q = ^ i exp(— iq • rj. It is important to point out 
that electron correlations are considered in the sys- 
tem through the many-particle dynamic structure fac- 
tor which is related to the density-density response func- 
tion X (q, w) through the fluctuation-dissipation theorem 
S(q,u)) = — (H/27rn s ) coth(hLu/2T)Im[x{q, u))], where n s 
is the electron density^ 

In the limit of low velocities u satisfying the condition 
Hq ■ u <C T, we may define an effective collision frequency 
as a proportionality factor between the momentum loss 
per unit time of the electron system and the average elec- 
ton velocity as 



dP 

= —mN e v\i. 

dt 



(23) 



The effective collision frequency is the sum of frequencies 
for a particular interaction Hamitonian v = y( e 9) -\-v( er ) + 
y{e%) _ rjijjg e i ec tron mobility is given then by 



j" = 



(24) 



Substituting the interaction Hamiltonians, given by 
Eqs. J5J) and ©, into Eq. J52J one obtains 



'.(eg) - 



-Xer) _ 



8m 3 T 



1 



AirmT 



q 3 S(q,0)dq; 



.Hqtanh(qd) , 
[ 2^ J 



(25) 



x\(xi\V rq (z)\xi}\ 2 N q S(q,to q )dq. (26) 

The calculation of is done in a similar way to that 
of T/( e s) and the result is straightforward 



-(ei) 



1 



AirmT 



2q 



x( X l\V lq (z)\ Xl )\ 2 S(q,0)dq. 



(27) 



In the BEFM approach, the form factor S(q, u>) is es- 
sential to the evaluation of Jl. Unfortunately S(q,uj) 
can be calculated analytically only in the case of the 
nonintcracting electron system and there are no reliable 
approximate expressions for S(q, to) in the whole fre- 
quency range in contrast to the static structure factor 
S(k) S(k, 0) which can be evaluated by appropriate 
approximation methods. In our calculations, we replace 
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S(q,u) in Eqs. Q25I271 by the noninteracting dynamical 
structure facto»2L2£ 



'27TTO 



1/2 



exp 



fi 2 g 2 mw 2 fiw\ 
8mT ~ 27V + 2T7 



and consider the limit Hu> q /T <C 1. 



(28) 



IV. RESULTS AND DISCUSSION 

Now we present the results of numerical calculations 
for SE mobility over a helium film from Eqs. |JSJ and (12411 . 
We have used the parameters £o an d I of the electron- 
interface scattering to have the best fit to the mobility 
experimental data. First we must say that the mobility 
curves in the BEA and FBEM exhibit the same overall 
behavior independently of the adjustable parameters for 
the electron-interface interaction. Results for the mo- 
bility dependence on the film thickness in the BEA are 
presented in Fig. 1 for a glass substrate (e s = 7.3) for 
T = 1.5 K and the fitting parameters £o = 10 A and 
I — 1420 A. For a sake of comparison, we plot the mo- 
bility data taken from Ref. [jJI One can see the deci- 
sive role of electron-interface scattering and the best fit 
is attainable for d < 300 A, where, as expected, the in- 
fluence of roughness is more pronounced. We also show 
separately the contributions to the mobility coming from 
electron-ripplon and electron-atom scattering processes. 
We observe that neither electron-gas nor electron-ripplon 
scattering mechanisms can explain the experimental data 
for wide range of d. By fitting the experimental data with 
mobility calculations in the FBEM, we obtain the best 
results for £ = 6.45 A and I = 5000 A. 

In Fig. 2 we depict, in a similar way to Fig. 1, mo- 
bility curves as a function of the film thickness for solid 
neon (e s = 1.19) and T = 1.2 K. The best fit of exper- 
imental dataiSi is achieved for £o = 180 A and I = 300 
A in BEA and £ = 85 A and I = 300 A in FBEM. We 
see now that £o and I are of the same order of magni- 
tude of d which makes unjustifiable the approximation 

(r) | jd <C 1 which supports the perturbative approach 
to obtain the potential Vi q {z) of Eq. J^J, and, hence, the 
applicability of Born approximation for the description 
of electron-interface scattering. We would guess that the 
neon surface, being much rougher than the glass one, it 
has more irregularities leading to significant fluctuations 
of d. Moreover, we observe a worse agreement between 
the theoretical curve and experimental data at large d, 
where the electron-ripplon and electron-gas mechanisms 
would dominated the scattering. 

The temperature dependence of SE mobility is de- 
picted in Fig. 3 where the same values for £o and I ob- 
tained by fitting /i(d) are used. As it can be seen we 
get rather good agreement with the experimental data of 
Ref. for d = 350 A. However the agreement becomes 
less satisfactory for smaller d. By considering only the 
electron-ripplon and electron-atom scattering we found a 
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FIG. 1: Surface electron mobility as a function of the film 
thickness for the helium film over a glass calculated within 
the BEA. The solid line represents the general result when all 
scattering processes are taken into account. The dashed line 
is the contribution from interface roughness scattering and 
the dotted line is the sum from contributions of the electron- 
ripplon and electron- gas scattering. The squares are experi- 
mental points of Ref. Ha 




FIG. 2: Same as in Fig. 1 but for a solid neon substrate. 
The solid line is the contribution coming from all scattering 
processes and dashed line represents the mobility when only 
the electron-interface roughness scattering is considered. The 
experimental points are taken from Ref. uu . 



difference of about more than one order of magnitude be- 
tween the calculation results and the experimental ones. 
Unfortunately, to our knowledge, there is no experimen- 
tal data on the temperature dependence of SE mobility 
in the case of a film over solid neon. 

Figure 4 shows the influence of the dielectric constant 
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FIG. 3: SE mobility as a function of temperature including 
all scattering processes over a glass substrate for some values 
of the film thickness. The marks correspond to experimental 
data taken from Ref. 

5 I 1 ' 1 ' 1 ' 1 1 1 ' 1 r- 



FIG. 5: Mobility as a function of the temperature calculated 
within the FBEM for d = 350 A and glass substrate. The 
solid line is the contribution from all scattering mechanisms 
whereas the dashed line is the contribution of gas scattering 
and the dotted line is the result for only the interface rough- 
ness scattering process. 




of the substrate on the SE mobility calculated within the 
BEA. We used here the same values of £o and / as in 
Fig. 1 for both glass and pmma (e s = 2.2) substrates. 
We see that the mobility varies inversely with respect to 
e s for whole of temperatures considered. This is a di- 
rect consequence from Ai ~ e s appearing in the collision 
frequencies which in turn are in the denominator of the 
mobility formulas [Eqs. (JSJ and Pijl]. 

In order to make clearer the role of different scatter- 
ing mechanisms we present in Fig. 5 the SE mobility, 
calculated within the FBEM, as a function of tempera- 
ture for a film over glass and d — 350 A. The calculated 
mobility curves are shown when the distinct scattering 
processes are considered separately. As it can be seen, 
the inclusion of roughness interface scattering is essential 
to explain the experimental data. The nonmonotonic de- 
pendence of SE mobility must be attributed to the cru- 



FIG. 6: Mobility as a function of the him thickness within 
the FBEM (solid line) and BEA (dashed line) for a glass sub- 
strate. 



cial role of electron-interface scattering since it is well 
known that the mobility limited only by electron-ripplon 
and electron-gas scattering is a monotonous decreasing 
function of temperature. 

Finally, we compare the mobility results obtained in 
the two different methods used here by depicting the 
curves of ^i(d) and /x(T) in Figs. 6 and 7 respectively 
for the case of a glass substrate. We plot also with the 
experimental data. Smaller values of £o obtained in 
FBEM as compared with those in BEA fulfill more satis- 
factorily the condition (r) | / <§; 1. We observe much 
better agreement between experimental data and theo- 
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FIG. 7: Mobility as a function of temperature. The notations 
are the same as in Fig. 6. 



rctical curves calculated in FBEM for the same values of 
£o and I. We point out that the numerical results within 
the FBEM are the same as those obtained in the com- 
plete control approximation (CCA)i22i2il In the CCA, the 
mobility is calculated with the BEA, but in the regime 
where the electron-electron collision frequency is much 
higher than the collision frequency due to other mecha- 
nisms. The electrons are supposed to redistribute their 
momenta, due to collisions between electrons, in such a 
way that the momentum of the total system does not 
change but electrons acquire the same drift velocity. The 
coincidence of results indicates that in order to establish 
the specific role of electron-electron interactions we must 



go beyond the noninteracting dynamical structure factor 
for calculating the mobility in the FBEM as given by Eq. 
Ij28jl . Our conclusions then are valid only in the low or 
intermediate electron density regime. 



V. CONCLUDING REMARKS 

In this work we have studied the dependence of the 
mobility of surface electrons over a liquid helium film on 
the temperature and on the film thickness. The trans- 
port properties are determined from the solution of the 
Boltzmann equation in the relaxation time approxima- 
tion and in the force balance equation. We employed 
the Gaussian correlated model to describe the interface 
roughness between the film and the substrate. The pa- 
rameters are adjustable through the evaluation of both 
temperature and film thickness dependencies of the mo- 
bility. The values obtained for the strength and range of 
the rough interface are of the order of few atomic layers 
and 100 A respectively which are reliable and in the meso- 
scopic regime. We showed that the SE mobility is limited 
by the roughness scattering because electron-ripplon and 
electron-gas atom scattering cannot explain the experi- 
mental data. However, we cannot rule out the presence 
of other mechanisms as SE localization by potential wells 
caused by the underlying solid substrate. 31 
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